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RESEARCH MEMORANDUM 


for the 

Bureau of Aeronautics^ Department of the Navy 
FLIGHT DETERMINATION OF THE LOW-LIFT DRAG AND LONGITUDINAL 

STABILITY OF A -i- - SCALE ROCKET-POWERED MODEL 
10 

OF THE DOUGLAS XF4D-1 AIRP.LANE AT 
MACH NUMBERS FROM 0.7 TO 1.4 
TED NO. NACA DE 34-9 

By Grady L. Mitcham^ Willard S. Blanchard^ Jr.^ 
and Earl C. Hastings^ Jr. 


SUMMARY 


A flight investigation has heen made to determine the drag and 
longitudinal stability of a i- scale model of the Douglas XF4D-1 air- 
plane from Mach numbers 0.7 to 1.4 at lift coefficients near zero. 

The drag rise occurred near M = 0.95* The external drag coeffi- 
cient was a constant value of about 0.012 at subsonic speeds up to the 
point of drag rise where it increased abruptly to a value of 0.030 at 
M = 1.0 followed by a more gradual increase to a value of O .038 at 
M = 1.25* The model indicated that^ at 35^000 feet and a level-flight 
free-stream Mach mrmber of 1.0, the drag of the full-scale airplane 
would exceed the thrust available from an XJ40-WE-8 engine with after- 
burning. The transonic trim change was small. The aerodynamic center 
moved gradually from the most forward location of 21.0-percent mean aero 
dynamic chord at M = 0.9 to the most rearward location of 40-percent 
mean aerodynamic chord at M = 1.25« The damping in pitch was low. 


INTRODUCTION 


An investigation of the drag and longitudinal trim characteristics 
of scale rocket -powered models of the Douglas XF4 d- 1 airplane at 
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-transonic and low supersonic speeds is being conducted by the Langley 
Pilotless Aircraft Research Division at the request of the Bureau of 
-Aeronautics^ Department of the Navy. 

The Douglas XF4D-1 is a jet-propelled^ single-place^ low-aspect- 
ratio^ svept-ving^ tailless^ interceptor -type airplane designed to fly 
at low supersonic speeds. The over-all investigation includes the 
effects of various combinations of external stores and rocket packages 
on the drag and longitudinal trim. 

The primary purpose of the test reported herein was to determine 
■' the drag at low lift coefficients for the airplane in the clean condi- 
tion. In addition to these data^ stability derivatives are also pre- 
sented since it was possible to analyze the short-period oscillations 
induced by booster separation and the trim change near a Mach mamber 
of 1.0. 


SYMBOLS 


M 

R 

V 

P 

H* 

q 

p 

A 

m 


S 

c 

W 


Mach number 
Reynolds number 
velocity^ feet per second 
static pressure^ pounds per square foot 
total pressure^ pounds per square foot 
dynamic pressure^ pounds per square foot 
density of air^ slugs per cubic foot 
duct area^ square feet 

mass flow through duct^ ^Pduct^duct“^duct^^ slugs per second 

mass of air flowing through a stream tube of area equal to 
the inlet-cowl area under free-stream conditions^ 
(Po^o^inlet slugs per second 

wing area of models square feet 

mean aerodynamic chord ^ feet 

weight of models pounds 
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longitudinal accelerometer reading 
normal accelerometer reading 

acceleration due to gravity^ 32.2 feet per second per second 


chord-force coefficient 


S Q,oS 


normal-force coefficient ( fiS. -H- 

\ S qoS 

drag coefficient (Drag/qQS) 
lift coefficient (Lift/qQS) 
thrust coefficient (Thrust/qQS) 


^hase ” 


hase-pressure coefficient of choking cup 


trim-normal -force coefficient 

rate of change of pitching-moment coefficient about center 
of gravity with respect to angle of attack, per degree 

moment of inertia in pitch, slug-feet^ 

period of short-period longitudinal oscillation, seconds 

time required for short-period longitudinal oscillation to 
damp to 1/2 amplitude, seconds 

rate of change of pitch angle with time, radians per second 

(i) 


rate of change of angle of attack with time, radians per 
second 
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Subscripts : 


0 

free 

inlet 

duct 

exit 

duct 

base 

base 


streson 

inlet 

exit 

of choking cup 


MODELS AND APPARATUS 
Models 


A three-view drawing of the scale model used in the present 

investigation is shown in figure 1. Photographs of the model are shown 
as figures 2 to 5* Table 1 gives the dimensional and mass characteris- 
tics. The model was constructed essentially of wood with aluminum 
inserts and castings . The eleven deflection was fixed ^ trailing edge 
up 0.3° prior to the flight test^ and the trimmer was built integral 
with the model and was not deflected. 

A cup designed to give choking at the exit was installed in the 
model. This installation was made in order that an accurate determina- 
tion of internal drag could be made with a minimum number of pressure 
measurements and to duplicate inlet Mach number of the full-scale air- 
plane at M = 1.4. A sketch of the choking-cup installation is shown in 
figure 1 and a photograph is shown as figure 4. The location of the 
total -pressure tube and the static-pressure orifice in one of the ducts 
is shown in figure 1. A plot of duct cross-sectional area from the inlet 
to the exit is shown in figure 6. 

Prior to flight testing^ the model was suspended by shock cords and 
shaken in the pitch plane with an electromagnetic shaker at frequencies 
up to 400 cycles per second. A fundamental frequency of 83 cycles per 
second (first bending) was observed from the telemeter record taken 
during the ground tests. As expected^ only the normal acceleration 
channel showed any frequency response in the ground tests. Resonance 
occurred at 83 ^ l 80 ^ and 240 cycles per second. 

The model was boosted to M = 1.4 by a solid -fuel_, 6-inch-diameter 
Deacon rocket motor which produces an average thrust of 650 O pounds for 
s-PP^^oximately 3*1 seconds. The model contained no internal -rocket - 
sustainer motor. 
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Launciiing vas accomplished from the zero -length launcher seen in 
figure 5* 


Apparatus 

The flight time history of the model was transmitted and recorded 
hy a telemeter system which gave eight continuous channels of informa- 
tion. The information measured during this flight test consisted of 
normal^ longitudinal^ and transverse acceleration^ free-stream total 
pressure^ inlet total pressure^ inlet static pressure^ exit static pres- 
sure^ and choking-cup hase pressure. A radiosonde released at time of 
firing was used to obtain free-stream temperature and static pressure. 
Ground equipment consisting of a CW Doppler radar unit and a radar 
tracking unit was used to determine model velocity and position in space 


METHOD OF ANALTCIS 


All data reported herein were obtained during the decelerating 
portion of the flight following separation of the model from the booster 

Drag . - Since the model flew at virtually zero normal-force coeffi- 
cient throughout the test speed range ^ the total-drag coefficient was 
assumed equal to the chord-force coefficient. The values of drag coeffi 
cient were ascertained by two independent methods. The first method 
involved the use of a longitudinal accelerometer mounted in the model 
and the second method involved differentiation with respect to time of 
the velocity along the flight path as determined from Doppler radar and 
tracking radar . 

The choking cup’ (figs. 1 and 4) was designed to maintain a Mach num 
ber of 0.32 at the duct inlets and a Mach number of 1.00 at the exit at 
supersonic speeds. With these conditions known^ mass flow and internal 
drag can be determined using exit static-pressure and free-stream condi- 
tions (see reference l) . Internal drag was determined by means of the 
relationship 

^^internal " ^ " '^exi^Pexit “ Po)^ 

Base drag was determined from 

_ -(Phase - Po)Abase 
Phase q^S 
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External drag was obtained by subtracting internal and base drag 
from total drag 

C ““ C -• C “ ^TV. 

^external ~ ^total ^internal ^ase 

Lift . - Since the lift coefficient throughout the flight was very 
low, lift was determined from the relationship 

Cl - % ^ 


where the term (an/g) was obtained from the telemeter signal trans- 
mitted by the normal accelerometer mounted in the model. 

Static longitudinal stability .- Oscillations in pitch were excited 
by several disturbances which occurred during the flight. The first 
oscillation was caused by the sudden change in trim when the model and 
booster separated^ the second was the result of the transonic trim 
change; the third was apparently caused by a gust which occurred at a 
Mach number of about 0.85 . These oscillations were analyzed for the 
determination of the period of the short-period longitudinal oscillation. 
Static longitudinal stability was then obtained as follows : 


■'“a “ " p 

pV S^ 


where was obtained from reference 2. 

^ Damping in pitch .- Where possible, the pitch oscillations were 

analyzed in an attempt to ascertain felt, however, that 

only the value obtained immediately after separation (highest Mach num- 
ber) is reliable quantitatively; the other two values obtained serve 
only as qualitative indications that damping was low, since the amplitudes 
of the oscillations were so small that accurate determination of 

1 . was impossible. Total damping was determined from the relationship 
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More complete discussions of the methods used in reducing data of this 
type are given in references 3 a^nd 4. 


PRECISION OF BkUk 


The extimated maximum errors in some of the data are presented in 
the following table: 



Supersonic 

Subsonic 

^total 

±0.0005 

±0,001 

internal 

+.0001 


^^ase 

±.0001 

±,0002 

^external 

±.0007 


m/mo 

±.0i<- 


Hinlet/Ho 

±.ok 

±.o8 

M 

±.005 

±,01 


RESULTS AND DISCUSSION 


The variation of Reynolds number with Mach number is shown in fig- 
inre 7<> The center of gravity of the model was located l6.7 percent 
behind the leading edge of the mean aerodynamic chord throughout this 
investigation . 


Trim Lift Coefficient 

The range of trim-lift coefficient of this flight investigation is 
shorn as a function of Mach number in figure 8. Since the elevens were 
fixed at a very small deflection (-0.3°)^ figure 8 represents the longi- 
tudinal trim characteristics of the airplane in the clean condition. 

The transonic trim change was small amounting to only ±0.025 change in 
lift coefficient between M = 0.95 and M = 1.02. 
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Buffet and Flutter 


As mentioned previously, resonant frequencies of the model were 
determined prior to the flight test. Analysis of the telemeter records 
of the flight, however. Indicated that no buffet or flutter was present 
during the flight. 


Drag 

The base-pressure coefficient of the choking cup is presented in 
figure 9. The results show that the base pressure decreases from a 
positive pressure at subsonic speeds to a negative pressure near M — 1.0 
and remains negative at supersonic speeds . This phenomenon has been 
observed on other models with comparable boattall angles at the base 
(reference 5) . This base pressure converted to base-drag coefficient is 
shown in figure 10. As might be expected, the base drag represents a 
small part of the total drag. 

The internal-drag coefficient (fig. 10) was roughly constant at 
0.0009 above M = 1.0. Values of internal -drag coefficient at subsonic 
speeds were not obtained since it was not possible to obtain a Mach num- 
ber of 1.0 at the duct exit. 

A summary of the drag data is given in figure 11. The total drag 
coefficient is nearly constant at a value of 0.012 from M = 0.72 to 
M = 0.95, the point at which the drag rise begins, followed by an abrupt 
increase to a value of 0.030 at M = 1.0, than a more gradual increase 
to a value of 0.039 at M = 1.25. Values of external-drag coefficient 
obtained from unpublished wind-tunnel data for the same configuration 
are shown plotted in figure 11. As indicated in the preceding discussion 
and figures, the contributions of base drag and internal drag to the 
over-all drag are small. These two Increments of drag were subtracted 
from the total drag coefficient to give the external-drag coefficient 
as presented in figure 11. Assuming that the internal-drag coefficient 
below M = 1.0 is equal to the internal-drag coefficient above M = 1.0, 
the values of total and external -drag coefficients are virtually the 
same at subsonic speeds . At supersonic speeds the external drag is 
smaller than the total drag by an increment in drag coefficient of about 
0.0012. Mass-flow ratios for this test and the unpublished wind-tunnel 
test are shown in figure 12. As might be expected, the higher mass-flow 
ratio corresponds to lower external drag. 

Also shown in figure 11 is estimated thrust coefficient of the 
XjiiO-'WE-S engine with afterburning at an altitude of 35^000 feet. The 
thrust coefficient intersects the zero lift-drag coefficient at M = 1.0^ 
therefore, when drag due to lift is considered the maximum level— flight 
Mach number will be less than 1.0. 
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Two factors possibly contributing to the high drag of this con- 
figuration are the sharp boattailing of the fuselage and the duct inlets 
The effect of comparable degrees of boattailing on the drag of bodies 
of revolution is shown in reference 6. Reference 7 indicates that this 
type of inlet tends toward a drag coefficient which increases with Mach 
number well into the supersonic range. 


Relative Mass Flow 

The variation of relative mass flow with Mach number is shown in 
figure 12. These values of relative mass flow were controlled by the 
choking cup located in the rear of the model. 


Total Pressure Recovery 

The total pressure in the duct was measured at station 24.30 as 
shown in figure 1. The measurements at this station were made to deter- 
mine the existence of instability of flow from two ducts discharging 
into a common duct. Mass flow based on the total -pressure and static- 
pressure tubes in the duct indicated no flow instability. The total- 
pressure recovery measured by the total-pressure tube at station 24.30 
is presented as a function of Mach number in figure 13. At the higher 
Mach numbers of the test the losses are considerably larger than those 
experienced through a normal shock. These losses are probably the 
result of flow separation in the vicinity of the duct inlet. 


Static Longitudinal Stability 

As stated in the introduction^ the primary purpose of this investi- 
gation was to determine drag at low lift coefficients. However^ as a 
result of several disturbances in pitch which occurred during the flight 
it was possible to measure the period of the short-period longitudinal 
oscillations. These values are presented in figure l4. The values for 
period were used to calculate the static longitudinal stability param- 
eter which is shown in figure 15. The center-of-gravity location 

of the model (l6.7 percent c) was forward of the center-of-gravity 
location quoted for the. full-scale airplane in reference 2. 

The values for slope of the lift curve given in reference 2 were 
used in the calculation of the aerodynamic -center location which is 
presented in figure l6. The aerodynamic center moves gradually from 
the most forward location of 21.0-percent mean aerodynamic chord at 
M = 0.9 to the most rearward location at 40-percent mean aerodynamic 
chord at M = 1.25. 
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Damping in Pitch 

The short-period longitudinal oscillations which occurred during 
the flight (see ’'Method of Analysis”) were analyzed to obtain a limited 
amount of damping data. The maximum magnitude of the oscillations 
decreased with decreasing Mach number^ thereby decreasing the accuracy 
of shown in figure 17. Values of total damping factor 

Cm^ + Cm^ based on these values of T 1/2 are shown in figure I 8 . 

Although these values are of limited accuracy, as seen in figure I 8 , 
they do indicate that damping in pitch at low lift coefficients is low 
and that further investigation is warranted. 


CONCLUSIONS 


From an analysis of the results of the flight test of a scale 

rocket -powered model of the Douglas XF4D-1 airplane from M = 0.7 to 
M = 1.4 the following conclusions are indicated: 

1. The external drag coefficient was approximately a constant value 
of 0.012 from M = 0.75 to M = 0.95^ the point at which the drag rise 
began, then increased abruptly to a value of 0.030 at M = 1.0 followed 
by a more gradual increase to a value of O.O 38 at M = 1.25- 

2. At 35,000 feet, maximum level-flight Mach number of the full- 
scale airplane will be less than 1.0 when the Xj40-WE-8 engine with 
afterburning is used. 

3 . The transonic trim change was small. 

4. The aerodynamic center moved gradually from the most forward 
location of 21.0-percent mean aerodynamic chord at M = 0.9 "to the 
most rearward location of 40-percent mean aerodynamic chord at M = 1.25. 
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5. The damping in pitch at lift coefficients near zero was low. 
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TABLE I 

PHYSICAL CHARACTERISTICS OP A ^ - SCALE MODEL OF XF4D-1 AIRPIANE 


Wing: 

Area (included), sq ft 5 *57 

Span, ft 3.35 

Aspect ratio 2.01 

Mean aerodynamic chord, ft 1.82 

Sweepback of leading edge, deg 52.5 

Dihedral (relative to mean thickness line), deg 0 

Taper ratio (Tip chord/Root chord) 0.33 


Airfoil section at root NACA 0007-63/30-9.5° Mod. 

Airfoil section at tip KACA 0004. 5-63/30-9 *5° Mod. 


Vertical tail: 

Area (leading edge extended to center line), sq ft 0.^ 

Aspect ratio . 2.o8 

Height (above fuselage center line), ft . 1.00 

Sweepback of leading edge, deg 33 

Taper ratio (Tip chord/Root chord) 0.26 

Airfoil section at root 1 . . . NACA OOO 8 - 63 / 3 O- 9 ® 


Airfoil section at tip NACA 0006-63/30-6.45° 


Elevon : 

Area (one), sq ft 0.23 

Span (one), ft 1.12 

Chord, ft 0.22 


Weight and balance: 

Weight, lb 109.9 

Wing loading, Ib/sq ft 19-7 

Center -of -gravity position (percent c) I6.7 

Moment of inertia in pitch, slug/ft^ 3*90 
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Figure 1.- Three-view drawing of the model; all dimensions are in inches 
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Figure 4 .- View showing installation of the choking cup. 
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Figure 5.- The booster-model combination on the mobile launcher. 
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Figure 6.- Duct areas. 
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Mach number, M 

Figure 7 *- Reynolds number. 
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Figure 9«- Pressure coefficient of the choking-cup base 
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Figure 10 . - Internal and "base drag. 
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Figure 11.- Lift, drag, and thrust summary. 
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c.g. at .167c 
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Figure l4.- Period of the short-period longitudinal oscillation. 
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Figure 15.- Rate of change of the pitching-moment coefficient with 

angle of attack. 
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Figure l6. - Aerodynamic center. 
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